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ABSTRACT 
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X 



Aims. In this paper we investigate the deceleration, jerk and snap parameters to distinguish between the dark energy and modified gravity 
models by using high redshift gamma-ray bursts (GRBs) and supernovae (SNe). 

Methods. We first derive the expressions of deceleration, jerk and snap parameters in dark energy and modified gravity models. In order to 
constrain the cosmographic parameters, we calibrate the GRB luminosity relations without assuming any cosmological models using SNe 
la. Then we constrain the models (including dark energy and modified gravity models) parameters using type la supemovae and gamma-ray 
bursts. Finally we calculate the cosmographic parameters. GRBs can extend the redshift - distance relation up to high redshifts, because they 
can be detected to high redshifts. 

Results. We find that the statefinder pair (r, s) could not be used to distinguish between some dark energy and modified gravity models, 
but these models could be differentiated by the snap parameter Using the model-independent constraints on cosmographic parameters, we 
conclude that the ACDM model is consistent with the current data. 
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CN 1- Introduction 



Recent observations of the Hubble relation of distant Type la 
supernovae (SNe la) have provided strong evidence for acceler- 
, ation of the present universe (Riess et al. 1998; Perlmutter et al. 
' 1999). The observations of the spectrum of cosmic microwave 
. background (CMB) anisotropies (Spergel et al. 2003;2007), 
' large-scale structure (LSS) (Tegmark et al. 2004; Eisenstein 
et al. 2005) and the distance-redshift relation to X-ray galaxy 
clusters (Allen et al. 2004; 2007) also confirm that the uni- 
verse is accelerating. Possible explanations for the accelera- 
tion have been proposed. A negative pressure term called dark 
energy is taken into account, such as the cosmological con- 
stant model with equation of state w - p/p = -1 (Weinberg 
1989), an evolving scalar field (Feeble & Ratra, 1988, Caldwell 
et al. 1998), the phantom energy for which the sum of the 
pressure and energy density is negative, and the Chaplygin 
gas (Kamenshchik et al. 2001). All the above models for ac- 
celeration are obtained by introducing a new energy compo- 
nent called dark energy. Alternative models, in which grav- 
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ity is modified, can also drive the universe acceleration, e.g., 
the Dvali-Gabadadze-Pon-ati (DGP) model (Dvali et al. 2000; 
Deffayet et al. 2002), Cardassian expansion model (Freese & 
Lewis 2002; Wang et al. 2003), and the f(R) gravity model 
(Vollick 2003; CaiToll et al. 2004). 

These two families of models, dark energy and modified 
gravity, are fundamentally different. An important question is 
whether it is possible to distinguish between the modified grav- 
ity and dark energy models that have nearly the same cosmic 
expansion history. Many works have been done on this topic. 
A usually-discussed quantity is the growth rate of cosmological 
density perturbations, which should be different in the models 
depending on different gravity theory even if they have an iden- 
tical cosmic expansion history. Recently, there have been ex- 
tensive discussions on discriminating dark energy and modified 
gravity models using the matter density perturbations growth 
factor (Linder 2005). But Kunz and Sapone (2007) demon- 
strated that the growth factor is not sufficient to distinguish 
between modified gravity and dark energy (Kunz & Sapone 
2007). They found that a generalized dark energy model can 
match the growth rate of the Dvali-Gabadadze-Porrati model 
and reproduce the 3-1-1 dimensional metric perturbations. 
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On the other hand, the statefinder pair (r, s) has also been 
proposed to distinguish between the models, where r = a/aH^ 
and s = (r - l)/3iq - 1/2). Sahni et al. (2003) demonstrated 
that the statefinder diagnostic could effectively discriminate 
different forms of dark energy (Sahni et al. 2003). Alam et al. 
(2003) investigated the cosmological constant, quintessence, 
Chaplygin gas, and braneworld models using the statefinder di- 
agnostic, and found that the statefinder pair could differentiate 
these models (Alam et al. 2003). Diff^erent cosmological mod- 
els exhibit qualitatively difl'erent trajectories of evolution in the 
r-s plane. The statefinder diagnostic has been extensively used 
in many models (Gorini et al. 2003). But the statefinder pair 
is difiicult to measure by cosmological observations (Visser 
2004; Cattoen & Visser 2007). The present values of cosmo- 
graphic parameters can be determined from observations (Riess 
et al. 2004; Visser 2004). Caldwell & Kamionkowski (2004) 
showed the jerk parameter could probe the spatial curvature 
of the universe (Caldwell & Kamionkowski 2004). The de- 
celeration, jerk and snap parameters are related to the sec- 
ond, third and fourth derivative of the scale factor respectively. 
Visser (2004) expanded the Hubble law to fourth order in red- 
shift including the snap parameter and put constraints on the 
deceleration and jerk parameters using SNe la (Visser 2004). 
Rapetti et al. (2007) constrained the deceleration and jerk pa- 
rameters from SNe la and X-ray cluster gas mass fraction mea- 
surements. For a redshift range of SNe la the terms beyond 
the cubic power of Hubble law can be neglected. In order to 
put a narrow constraint on the snap parameter, we need high- 
redshift objects. GRBs may be a useful tool. GRBs can be 
detectable out to very high redshifts (Ciardi & Loeb 2000). 
The farthest burst detected so far is GRB 090423, which is at 
z = 8.2 (Olivares et al. 2009). A lot of work in this so-called 
GRB cosmology has been pubhshed (Dai, Liang & Xu 2004; 
Ghirlanda et al. 2004; Di Girolamo et al. 2005; Firmani et al. 
2005; Friedman & Bloom 2005; Lamb et al. 2005; Liang & 
Zhang 2005, 2006; Xu, Dai & Liang 2005; Wang & Dai 2006; 
Schaefer 2007; Wright 2007; Wang, Dai & Zhu 2007; Gong & 
Chen 2007; Li et al. 2008; Liang et al. 2008; Qi, Wang & Lu 
2008a,b; Basilakos & Perivolaropoulos 2008; Kodama et al. 
2008; Wang, Dai & Qi 2009). Very recently, Schaefer (2007) 
used 69 GRBs and five relations to build the Hubble diagram 
out to z = 6.60 and discussed the properties of dark energy in 
several dark energy models (Schaefer 2007). He found that the 
GRB Hubble diagram is consistent with the concordance cos- 
mology. Liang et al.(2008) calibrated the luminosity relations 
of GRBs by interpolating from the Hubble diagram of SNe la 
at z < 1.4 with the assumption that objects at the same redshift 
should have the same luminosity distance (Liang et al. 2008). 
This method is model-independent. More recently, Capozziello 
& Izzo (2008) used the Liang et al. (2008) results to constrain 
the cosmographic parameters and found the results cahbrated 
by SNe la data, agree with the ACDM model. Cardone et al. 
(2009) used 83 GRBs and six correlations to build the Hubble 
diagram. Butler et la. (2009) found a real, intrinsic correlation 
between Ei^o and fipeak using latest Swift GRB sample. 

Riess et al. (2004) found that the jerk jo is positive at the 
92% confidence level based on their "gold" dataset and is pos- 
itive at the 95% confidence level based on their "gold-i- silver" 



dataset. Neither explicit upper bounds are given for the jerk nor 
are any constraints placed on the snap s^- Rapetti et al. (2007) 
measured = -0.81 ± 0.14 and j = 2A6ll f^ ^ ^ model 
with constant jerk (Rapetti et al. 2007). Capozziello & Izzo 
(2008) used 27 GRBs to derive the values of the cosmographic 
parameters. They found = -0.78 ± 0.20, 70 = 0.62 ± 0.86 
and So = 8.32 ± 12.16. In this paper, we use more GRB data to 
constrain the cosmography parameters in several dark energy 
and modified gravity models. 

In this paper, we calibrate the luminosity relations of GRBs 
using SNe la and calculate the deceleration, jerk and snap pa- 
rameters of several dark energy and modified gravity models 
using SNe la and GRBs. We also use a model-independent 
method to constrain the cosmographic parameters. We find that 
in some models the jerk parameter is almost equal to each other. 
So this parameter is not used to distinguish between the mod- 
els. However, the snap parameter in all the models is different, 
so we can distinguish between the models using the snap pa- 
rameter. 

The structure of this paper is organized as follows. In sec- 
tion 2 we introduce the bubble, deceleration, jerk and snap pa- 
rameters. In section 3 we derive expressions of cosmographic 
parameters of the Hubble law in several dark energy models. 
In section 4 we present expressions of cosmographic parame- 
ters of the Hubble law in modified gravity models. The con- 
straints on models parameters and cosmographic parameters of 
the Hubble law are given in section 5. Finally, section 6 con- 
tains conclusions and discussions. 

2. Hubble, deceleration, jerk and snap parameters 

The expansion rate of the Universe can be written in terms of 
the Hubble parameter, H = a/ a, where a is the scale factor and 
d is its first derivative with respect to time. As we known that 
q is the deceleration parameter, related to the second derivative 
of the scale factor, j is the so-called "jerk" or statefinder pa- 
rameter, related to the third derivative of the scale factor, and s 
is the so-caUed "snap" parameter, which is related to the fourth 
derivative of the scale factor. These quantities are defined as 
1 a 



■- 

- lA 



(1) 



(2) 



(3) 



The deceleration, jerk and snap parameters are dimensionless, 
and a Taylor expansion of the scale factor around to provides 

ait) = ao |l + Hoit - to) - ^qoHlit - tof + ^ joH^it - to)^ 

+ ^soH'o(t-to)UO[{t-tof]K4) 
and so the luminosity distance 

= 1^ + ^(1 - qo)z^ -^[l-qo-3ql + jo)z^ 
+^ [2 - 2^0 - 15^^ - 15^^ + + 10^070 + so] + 0(z^))(5) 
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(Visser 2004). For the redshift range of SNe la the terms be- d£ 

yond the cubic power in Eq. (5) can be neglected. If models dz o 

have the same deceleration and jerk parameters, we can see ^ 

degeneracy of these models from Eq. (5). Therefore we must j^'^^ = -(2 + 9(1 - Qm)wo + 9(1 - Qm)Wo)' 
measure the snap parameters to distinguish between the mod- 
els. This needs high-redshift objects. The relations among the Y^p,^^ 1 
q(z),jiz) and ^(z) are 
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(16) 
(17) 



_XCDM 
*0 



= _(_14 - 81(1 - Qm)wo - 9(16 - 19Qm + ^Q^iH 



j(z) = q(z)+2q\z) + (l+zAzy, 

dz 

s(z) - -(1 + z)^(z) - 2j(z) - 3j(z)q(z). 
dz 

The Friedmann equation is 



(6) 

(7) 

(8) 



From Einstein's equations, we can obtain the dynamical equa- 
tion of universe 



a AnG v-i , ^ „ , 

i 

The conservation equation is 
Pi + 3H(pi + Pd = 0. 



(9) 



(10) 



a = — — > [dpi(l + 3wj) + api(l + 3w,) + ap, X 3w,] 



AnG 

= — 2^[aHpiil + 3wi) - 3Hapiil + w,)(l + 3w,) 



+api X 3^,], 



(11) 



AnG v-i d[dpi(l + 3wj) + ap,(l + 3w,) + ap,- X 3^,] 
^ = -—l^ dt 

AnG ^ 

= — — 2_Jdpi{l + 3wi) + 2dpi{\ + 3wi) + 6apiW,- 
-i-ap,(l + 3w,) -I- 6ap,Wj -i- api x 3w,]. 



(12) 



3. Dark energy models 

3.1. w{z) parameterization model 



w(z) = Wo 

For this model, we obtain 



2[l+^o(l-i^M)]-l, 



(14) 



dq 

Tz 



XCDM 



-27(3-4tlM + n^)w^). 



(18) 



These expressions are consistent with Bertolami & Silva 
(2006). 

A more interesting approach to explore dark energy is to 
use time-dependent dark energy model. The simplest parame- 
terization including two parameters is (Maor et al. 2001 ; WeUer 
& Albrecht 2001) 



w{z) -Wq + WiZ. 



(19) 



In this dark energy model the luminosity distance is (Linder 
2003) 



dL = cHq\1 +z) 



f 

Jo 



dz[(l+Z?nM 



\3(l+M'o-»'i)„3wiz-|-l/2 



^WCDM 
^0 



I 3 

- + -(1 - ^m)wo 



(20) 
(21) 



In order to derive the jerk and snap parameter, we differentiate ■'o 
Eq.(9) 



.•WCDM 



1 



- (2 + 9( 1 - Qm) WO + 9( 1 - £Im)wI + 3(1- Qm)wi )(22) 



„WCDM 
*0 



1 



(-14 - 9(16 - 19Qm + 3Qt)w, 



-27(3 - AQ.M + ^m)wI - 45(1 - Qm)wi 
+9wo(l - Qm)(-9 - (Qm - 7)wi)) 



(23) 



We consider the ChevaUier-Polarski-Linder parameteriza- 
tion (Chevallier & Polarski 2001; Linder 2003) 



w(z) = Wo -I- 



l+z' 



(24) 



The luminosity distance is (Chevallier & Polarski 2001; Linder 
2003) 

dL = cHq\\+z} r dz[{l+z?nM 
Jo 

+(1 - Qm)(1 + z)3(l+"'o+wi)g-3wiz/(l+dj-l/2_ 

The cosmographic parameters are: 



(25) 



1 3 

2 2*-^ 



(26) 



We first consider the dark energy with a constant equation of jCpl 
state. 



1 



-(2 + 9(1 - Qm)wo + 9(1 - ^m)wI + 3(1 - Qm)wi), (27) 



(13) rpL _ 1 



-(-14-9(16-19Qm + 3Q1,)vw^ 

-27(3 - A€Im + O^M>l - 33(1 - ^m)wi 
+9wo(l - Om)(-9 - (Qm - 7)wi)). 



(28) 



Capozziello, Cardone & Salzano (2008) and Capozziello & 
Izzo (2008) also derived cosmographic parameters in this 
^^^^ model. Our results are consistent with theirs. 
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3.2. generalized Chaplygin gas model 

We consider the generalized Chaplygin gas (GCG) model, 
which is characterized by the equation of state 



equation is modified as Eq. (39). The universe undergoes ac- 
celeration requires n < 2/3. If n = 0, it is the same as the 
cosmological constant universe. We can obtain H(z) by using 
Eq. (O andp,„ = p,„(l + = ^„Pc(l + z)\ 



PGCG = -A/PccG- 



(29) 



HizY = m[^l„,(l + zY + (1 - "mXl + zf"]. 



(40) 



We can integrate the conservation equation for generalized 
Chaplygin gas, leading to 



PGCG = Pgcgo[As + (1 - A,)a 



3(l+Q-)-|l/(l+a) 



(30) 



where pcho is the energy density of GCG today, and Aj = 
A/p^^y The attractive feature of the model is that it can unify 
dark energy and dark matter. The reason is that, from Eq. (30), 
the GCG behaves as dustlike matter at an early epoch and as 
a cosmological constant at a later epoch (Kamenshchik et al. 
2001; Bento et al. 2002). The Friedmann equation can be ex- 
pressed as 



where pc - 3Hq/^7:G is the critical density of the universe. The 
luminosity distance in this model is 

dL^cHo\l+z) I dz[(l+zM„ + (l-^l,„)(l+zY"]-"^.(4l) 



For the Cardassian expansion model, we obtain 

1 3 
- + -( 

2 2 



^ - + ^(l - n)(^M - I), 



dq 
Tz 



Card Q 

= -(«- 1)2(1 -Qm)Qm, 
^ 



H\z,Ho,A„a) = HlE\z,A„a), 
where 



E\z,A„a) = ni,(l+zy+(l-nh)[AMl-AMl+zf'^"^]^,02) 



(31) ^ 

dz 



di 27 

^—(n-lf(l-QM)^Mn, 
^ 



C^>d = -(2 + 9n(QM - 1) + 9n\l - Qm)) 



(42) 
(43) 
(44) 
(45) 



Qi is the density parameter of the baryonic matter. The lumi- 
nosity distance is 



Jo 



di = cH-'{l + z) dz{{l + zY^h + (1 - ^b) 



[A, + (l-A,)(l+z)3(i+'"] 



-1/2 



„Card 



= ^(4- 18QM-27n^(3-4QM + Ol,) 



-9n(4 - 1Q.M + 30i,) + 9n\l 1 - 17Qm + 6n^))(46) 



(33) 4.2. Dvali-Gabadadze-Porrati model 



For the GCG model we obtain (Bertolami & Silva 2006; Wang, 
Dai & Qi 2009) 



^GCG 



= :^{l-As)-\ 



dz 

</«CG 

dz 



= -A,(l-A,)(l+a), 

^ 



= -ya(l+a)(2A,-l)(A,-l)A, 



Jo''° = 4(1-^.^X1+ (3 +6t^)A,), 



(34) 
(35) 
(36) 
(37) 



In the DGP model the modified Friedmann equation due to the 
presence of an infinite-volume extra dimension is (Deffayet et 
al. 2002) 



where the bulk-induced term, Qr^ , is defined as 



(47) 



(48) 



„GCG 



= -(A,-l)(7-H6(2 + a-6a'2)A,+9(-3+2ff-H8a2)A?).(38) 



For a flat universe, Q.^ - 0. In the above equation, is the 
crossover scale beyond which the gravitational force follows 
the 5-dimensional Ijr' behavior Note that on short length 
scales r «; (at early times) the gravitational force follows 
the usual four-dimensional l/r^ behavior. For a spatially flat 
universe, Q,. = (1 - Q.,y,flA. We obtain 



4. Modified gravity models 

4.1. Cardassian expansion model 

The original Cardassian model was first introduced in (Freese 
& Lewis 2002) as a possible alternative to explain the accel- 
eration of the universe. They modified the Friedmann equation 
as 



DGP _ i , ^ - 1 

^» "2 21+Qm^ 



dq 
Tz 

dj_ 

dz 



H^^'-fp,„^Bpl 



(39) 



This model has no energy component besides ordinary matter 
If we consider a spatially flat FRW universe, the Friedmann 



;DGP _ 
JQ 



DGP 



9Qm(1-^^m) 
(I+^mY ' 

^ 54Q^(1 - Qm) 
(1 + ^mY ' 

1 + 3Qm 



6"^ 



{I+QmY 

1 - 4Qm - 35QL - 26Qi + 32Ql - SOnl 



M 



(1+Qm)5 



(49) 



(50) 



(51) 



(52) 



(53) 
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4.3. f{R) gravity 

f{R) gravity models, in which the gravitational Lagrangian is a 
function of the curvature scalar R, also can explain current cos- 
mic acceleration (Volhck 2003; Carroll et al. 2004; Capozziello 
et al. 2009). Poplawski (2006) derived a quite complicated ex- 



GRB can be calculated. So the luminosity distances and dis- 
tance modulus can be obtained. After obtaining the distance 
modulus of each burst using one of these relations, we use the 
same method as Schaefer (2007) to calculate the real distance 
modulus. 



pression of jerk parameter in /(/?) =R-^ (Poplawski 2006): = (^ /i,/ct^ )/(^ o-'f), 



(58) 



/ = [(0/ - 3/)(2/'2 + 0/'/" - 6//")(3O0YV2 + 100V"/'" 

2ffif"2 11 A.2 fi?: fii ii^A.2ffi2 fill Q^fiA .nAffi^ where the summation runs from 1 - 5 over the relations with 

available data, //, is the best estimated distance modulus from 

the i-th relation, and cTm. is the corresponding uncertainty. The 

(2</.V"* - ' -■ ^ ' 



+ 1890//"' + 1890/V'/'" + 1920//'V" - 162/V"' 



-261 fm - i2cf>^f'^ + m'f'f" - 750V'V/" - 120//' 
+ m<pff'f" - 162/V") X (30Y/"2 - 150//"2 - 8/'3 
+18/V" + 27//7" - 0/'V" + 0'/"/"' - 90//'/'" 
+18/V'")] X [(0/' - 3/)3(2/'2 + 0/'/" - 6ff"frK 

The snap parameter in this model is (Poplawski 2007) 



6/(0/ -2/) . 80/ - 15/ 

(2/'2 + 0/'/"-6//") ^ 4>f'-3f ' 



(55) 



ju = 51og(JL) + 25. 



The likelihood functions can be determined fmmx^ statistic. 



2 

XsNt 



+ (rt 



1/2 



uncertainty of the distance modulus for each burst is 

(54) 
The;i'^ value is 



(59) 



Poplawski (2007) calculated go = -0.67+°°^, jo = i-Ol^ofi 
and So = -0.22;|;q j^. These expressions of the jerk and snap 
parameters are only valid in Palatini variational principle. A 
generic formulae of cosmographic parameter are derived by 
Capozziello, Cardone & Salzano (2008) (for more details, see 
Equations (23)-(33) in their paper). They also gave the best 
fitted value: go = -0.55 + 0.38, jo = 1.0 + 5.4 and sq = 
-0.35 + 28.1 using SNe la. In this paper, we only use Poplawski 
(2006) as an example for the /(/?) gravity. 

5. Constraints from SNe la and GRBs 

Davis et al. (2007) fitted the SNe la dataset that include 60 
ESSENCE SNe la (WoodVasey et al. 2007), 57 SNe la from 
Super-Nova Legacy Survey (SNLS) (Astier et al. 2006), 45 
nearby SNe la and 30 SNe la detected by HST (Riess et al. 
2007) with the MCLS2K2 method. With the luminosity dis- 
tance (//, in units of megaparsecs, the predicted distance modu- 
lus is 



(56) 



where o-y is the dispersion in the supernova redshift (trans- 
formed to distance modulus) due to peculiar velocities, juo,, is 
the observational distance modulus, and cr^„ . is the uncertainty 
in the individual distance moduli. The confidence regions can 
be found through marginalizing the likelihood functions over 
Ho (i.e., integrating the probability density p oc exp~^ ''^ for all 
values of Ho). 

We use the calibration results obtained by using the inter- 
polation methods directly from SNe la data (Liang et al. 2008). 
The calibrated luminosity relations are completely cosmology 
independent. We assume these relations do not evolve with red- 
shift and are valid in z > 1.40. The luminosity or energy of 



2 

Xgrb 



N 

!=1 



(60) 



where yUflt,, and (r^^y are the fitted distance modulus and its er- 
ror. 

We combine SNe la and GRBs by multiplying the likeU- 
hood functions. The total yjjyg is^rfotai = AtInc +ArGRB- 
best fitted value is obtained by minimizing 

5.1. Constraints on cosmograpliic pararrieters 

In our analysis, we consider the flat cosmology. We use h = 
0.72 + 0.08 from the Hubble Space Telescope key projects 
(Freedman et al. 2001). Riess et al. (2009) used the old dis- 
tance ladder and observed Cepheids in the near-infrared where 
they are less sensitive to dust and found h = 0.742 ± 0.036. 

Let us first consider observational constraints on dark en- 
ergy models. In Fig.l, we show the distribution probabihty 
a function of D.m in the flat ACDM model from SNe la and 
GRBs. From this figure, we have Qm = 0.27 ± 0.04. The cos- 
mographic parameters in ACDM model are q'o = -1 + jOm, 
jo = 1.0 and so = 1 - |^^m- We can obtain ^o = -0.60 ± 0.06, 
70 = 1.0 and so = -0.22 ± 0.18. 

In Fig. 2 we present constraints on Q.m and w from Icr to 
3cr using 192 SNe la and 69 GRBs in the w = wq model. We 
measure Qm - 0.29;^!! !! and wq - -1.04'^?^ l^. The cosmo- 



(57) graphic parameters in the w = wo model are go - -0.61;|;q 



+0.38 
,60' 



70= 1.13t}:iOandso = -0.08:i;^5 

In Fig. 3 we present constraints on wo and wi from Icr to 3cr 
using 192 SNe la and 69 GRBs in the w = wq + wiz model. The 
values of parameters are wq = -1.14 ± 0.19 and wi = 0.63 ± 
0.44. The cosmographic parameters are go = —0.75 ± 0.21, 
70 = 2.21 ± 0.93 and so = -12.25 + 9.18. 

In Fig.4 we present constraints on wo and wi from Icr to 
3o- using 192 SNe la and 69 GRBs in the w = wo + wiz/(l + z) 
model. We measure wq = -1.22 
cosmographic parameters are ^o = -0.90 ± 0.33, 70 
and so = -25.52+g-33. 

Fig. 5 shows constraints on As and a from Icr to 3cr using 
SNe la and GRBs in the GCG model. The parameters are As = 



0.30 and wi = 1.6!{ 20. The 

-'■^-'-2.09 
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Fig. 1. Luminosity distance - redshift diagram. The circles are 
the GRBs. The soUd line is the result of our fitting 




Fig. 3. The same as Fig.2 but in the w = wo+wiz model. 





Fig. 2. Constraints on Q.m and w from Icr to 3cr using 192 SNe 
la in the w = wq model. 



Fig. 4. The same as Fig.2 but in the w = wq +wiz/(1 +z) model. 



O.lS'^Q^y The cosmographic parameters 



0.79 + 0.12 and a 

are qo = -0.695 ±0.18, 70 = 1-18!^ ^ and sq = -0.37;!: J;^^. 

In Fig. 6 we present constraints on SIm and n from Icr to 3(T 
using 192 SNe la and 69 GRBs in the Cardassian expansion 
model. We measure Om = 0.29 ± 0.1 1 and n = -0.07:;:^-^4. The 
cosmographic parameters are qo = -0.6T^-^q, jo = 135'^\--^^ 
and so = 0.361|^i. 

Fig.7 shows constraints on D.m using SNe la and GRBs in 
the DGP model. The value of Qm is = 0.20 + 0.02. The 
cosmographic parameters are qo - -0.50 + 0.04, jo - 0.83 + 
0.02 and So = -0.56 ±0.12. 

We directly use Eq.(5) to constrain the cosmographic pa- 
rameters. This analysis uses the FRW metric only, so we have 
not specified any gravitational theory yet. The luminosity dis- 



tance only depends on redshift z and cosmographic parame- 
ters. So this method is fully model independent. We use the 
192 SNe la and 69 GRBs and find the best fit parameters are 
^0 = -0.85 + 0.19, jo = 1.50 + 0.80 and so = 3.49 + 4.50. The 
results are consistent with the flat ACDM model. 

In Table 1 we summarize the constraints on cosmographic 
parameters. The deceleration and jerk parameters in the w = 
Wo, GCG, Cardassian expansion and f(R) models are almost 
the same in the Icr confidence level. These values are consistent 
with the deceleration and jerk parameters of the ACDM model 
in the Icr confidence level. So these models can not be discrim- 
inated using the present value of the statefinder pair. However 
the snap parameter in all the models is difi'erent and thus can 
be used to discriminate the cosmological models. In the future. 
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Fig. 7. The same as Fig.l but in the DGP model. 



Fig. 5. The same as Fig.2 but in the GCG model. 




0.0 0.2 



0.4 0.6 



and gamma-ray bursts. GRBs can extend the redshift - distance 
relation up to high redshifts. We find that the deceleration and 
jerk parameters in the w - vvo, GCG, Cardassian expansion and 
f(R) models are almost the same in the Icr confidence level. So 
these models can not be discriminated using the present value 
of the statefinder pair. We find that the dark energy models and 
modified gravity models could be distinguished between by the 
snap parameter. Using the model-independent constraints on 
cosmographic parameters, we find the ACDM model is consis- 
tent with the current data. 
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Fig. 6. The same as Fig.2 but in the Cardassian expansion 
model. 

more data will give a precise snap parameter in difi'erent mod- 
els. 

6. conclusions and discussions 

The cosmic acceleration could be due to a mysterious dark en- 
ergy, or a modification of general relativity (modified gravity). 
In this paper we investigate the deceleration, jerk and snap pa- 
rameters in modified gravity models and dark energy models. 
We caUbrate the GRB luminosity relations without assuming 
any cosmological models using SNe la. Because gamma-ray 
bursts can be detected in high redshifts, we calculate the de- 
celeration, jerk and snap parameters using type la supemovae 
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model ^0 Jo 

ACDM -0.60 ±0.06 1.0 ±0.11 -0.22 ±0.18 



w 



-'-0.60 '-■'-•'-{.19 "-""-a.TS 

w = vvo + wiz -0.75 ±0.21 2.21 ±0.93 -12.25 ±9.18 

m; = Wo + -0.90 ± 0.33 3.93!^-^^ -25.52+g-33 

GCG " -0.70 ±0.18 1.18:°-^^ -0.37![-^'5 

Cardassian -Q-€n*^'^l ^2Ql\f^ 0.22:2|i 

DGP -0.50 ±6.04 0.83 ±6.02 -0.56 ± 6.12 

f(R) -0.67!°:°^ i.or»:°j -o-22!iii^ 



